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Abstract 

Evolution has contributed to the successful invasion of exotic plant species in 
their introduced ranges, but how evolution affects particular control strategies is 
still under evaluation. For instance, classical biological control, a common strat- 
egy involving the utilization of highly specific natural enemies to control exotic 
pests, may be negatively affected by host hybridization because of shifts in plant 
traits, such as root allocation or chemical constituents. We investigated introgres- 
sion between two parent species of the invasive shrub tamarisk (Tamarix spp.) in 
the western United States, and how differences in plant traits affect interactions 
with a biological control agent. Introgression varied strongly with latitude of 
origin and was highly correlated with plant performance. Increased levels of 
T. ramosissima introgression resulted in both higher investment in roots and 
tolerance to defoliation and less resistance to insect attack. Because tamarisk 
hybridization occurs predictably on the western U.S. landscape, managers may be 
able to exploit this information to maximize control efforts. Genetic differentia- 
tion in plant traits in this system underpins the importance of plant hybridization 
and may explain why some biological control releases are more successful than 
others. 



Introduction 

Rapid evolution is an important process contributing to 
the success of some invasive plant species in their intro- 
duced ranges (for reviews, see Lee 2002; Bossdorf et al. 
2005; Keller and Taylor 2008; Prentis et al. 2008). Some 
introduced plant species have evolved faster growth, higher 
fecundity, or altered allocation of herbivore defenses com- 
pared to populations in the native ranges (Blossey and 
Notzold 1995; Hull-Sanders et al. 2007). Moreover, evolu- 
tion is invoked to explain why populations of several exotic 
plant species persist at low numbers for decades before 
undergoing exponential population growth (e.g., termed 
'evolution of invasiveness', Ellstrand and Schierenbeck 
2000). Invasive species provide a unique opportunity to 
study evolution in contemporary time; it is also imperative 
to understand the evolutionary dynamics that influence 
invasion success because these organisms are related to 
declines in biodiversity, drastic changes in ecosystem func- 



tion, and costs of over $120 billion annually in the United 
States alone (Wilcove et al. 1998; Mack et al. 2000; Pimen- 
tel et al. 2005). 

Evolutionary mechanisms, such as founder events, intra- 
and interspecific hybridization, and adaptation, have all 
resulted in significant evolutionary change in the intro- 
duced ranges (Lee 2002; Bossdorf et al. 2005; Prentis et al. 
2008). For instance, despite several deleterious effects, 
hybridization can also result in novel genotypes which are 
better suited to their environment than either parent spe- 
cies (Lee 2002; Donovan et al. 2010). Hybridization has 
been implicated in the transfer of beneficial genes for traits 
such as cold hardiness or resistance to fungal diseases and 
herbivores (Snow et al. 1999; Abbott et al. 2003; Whitney 
et al. 2006; Rieseberg et al. 2007). Selection can act on 
these novel hybrids, spreading beneficial alleles rapidly 
throughout populations of invasive plants (Keller and Tay- 
lor 2008). Indeed, hybridization has been implicated in 
numerous plant invasions (Ellstrand and Schierenbeck 
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2000; Schierenbeck and Ellstrand 2009), with several hybrid 
taxa evolving larger size or higher fecundity than either 
parent species or outcompeting and replacing parent spe- 
cies (Campbell et al. 2006; Whitney et al. 2006; Ridley and 
Ellstrand 2010). Invasive plant hybrids indirectly threaten 
native communities and pose difficulty to land managers 
who are responsible for their control (Vila et al. 2000; Blair 
et al. 2008). 

One effective tool for managing invasive plants is classi- 
cal biological control, where specialized natural enemies 
are imported and released to provide top-down effects on 
their hosts (DeBach 1964; van Klinken and Raghu 2006). 
How biological control is affected by invasive plant hybrid- 
ization and other evolutionary mechanisms remains largely 
unknown (Muller-Scharer et al. 2004). In native systems, 
herbivores can distinguish among hybrid genotypes (Fritz 
et al. 1998; McGuire and Johnson 2006), and hybrid plants 
can have more or less resistance to attack compared with 
their parent species (Whitham 1989; Fritz et al. 1999; Krebs 
et al. 2011). In invasive plant systems, hybridization can 
affect the frequency of herbivore attack, especially when 
hybrids are compared as a group to their parent species 
(Blair et al. 2008; Krebs et al. 2011; Cuda et al. 2012). 
However, little is known about how susceptibility to her- 
bivory varies across levels of species introgression in inva- 
sive plants. This information could be useful in cases where 
invasions are comprised almost entirely of hybrid geno- 
types (Williams et al. 2007; Gaskin and Kazmer 2009). 
Clearly, more research is needed to determine the extent to 
which invasive plant hybridization affects the efficacy of 
classical biological control. 

In addition to hybridization, another mechanism of 
evolution important to invasive plants is adaptation to 
abiotic stress (Bossdorf et al. 2005; Ridley and Ellstrand 
2010). Introduced plants are exposed to novel climatic 
and edaphic conditions in the new range on which selec- 
tion can act. One way this selection is evident is through 
the evolution of latitudinal clines. When grown in a com- 
mon environment, populations of several exotic plant spe- 
cies exhibit inherited genetic differences in traits, such as 
size, phenology, and cold hardiness, reflecting the climate 
from the latitude where they were collected (Weber and 
Schmid 1998; Kollmann and Banuelos 2004; Maron et al. 
2004; Leger and Rice 2007; Montague et al. 2008; Keller 
et al. 2009). While latitudinal clines of native plant species 
have been well documented for quite some time (Turreson 
1930; Clausen et al. 1940), similar patterns are also seen 
for non-native plants, some of which have evolutionary 
histories as little as 100 years in their new ranges (Ridley 
and Ellstrand 2010; Hodgins and Rieseberg 2011). It is 
unclear whether evolution of strategies to cope with abi- 
otic stress may interact with biological control. For 
instance, common garden experiments by Friedman et al. 



(2008, 2011) show that northern populations of Tamarix 
spp. have evolved increased cold hardiness due to extreme 
minimum temperatures and have likely adapted to winter 
dieback by allocating more resources to belowground tis- 
sues (e.g., coarse roots). In turn, we hypothesize that these 
same individuals may have increased tolerance to leaf her- 
bivory because they have more carbohydrate reserves 
available for leaf flush following defoliation. In short, evo- 
lution of life history traits to cope with abiotic stress could 
also influence ecological interactions with specialized nat- 
ural enemies, consequently affecting the success of biologi- 
cal control efforts. 

Here, we investigate the effects of hybridization on plant 
traits and latitudinal clines of the perennial shrub tamarisk 
(i.e., saltcedar, Tamarix spp., family Tamaricaceae). Several 
species of tamarisk were introduced to the United States in 
the nineteenth century to stabilize stream banks and to 
serve as ornamental plants (Robinson 1965). By the mid- 
1900s, at least four species were considered serious pests in 
the arid west (Gaskin and Schaal 2002; Gaskin et al. 2012). 
Tamarisk has been estimated to occupy at least 360 000 ha 
in the western United States (Nagler et al. 2011) and is the 
second most dominant woody riparian species (in terms of 
percent cover) in the western United States (Friedman 
et al. 2005). In 2001, the biological control agent Dio- 
rhabda carinulata Brulle (Coleoptera: Chrysomelidae) was 
released at field sites across western North America. In 
some locations, populations of this insect became well 
established and control of tamarisk is being achieved, while 
at other sites, the control agent has failed to establish 
despite repeated attempts (DeLoach et al. 2003, 2009). 
Some authors have speculated that genetic makeup of host 
plants may play a role in these instances (Gaskin and Schaal 
2002; Gaskin and Kazmer 2009; Dudley et al. 2012; Hultine 
et al. 2013). Could differences in susceptibility among 
hybrid Tamarix genotypes be driving the success or failure 
of tamarisk biological control? 

A hybrid 'swarm' of Fl, F2, and backcrosses to two par- 
ent species, T. ramosissima and T. chinensis, make up the 
bulk of the tamarisk invasion in the western United States 
(Friedman et al. 2008; Gaskin and Kazmer 2009). Based on 
diagnostic markers and AFLP data, Gaskin and Kazmer 
(2009) estimated 83-87% of genotypes collected across sev- 
eral sites in western North America were indeed T. ramo- 
sissima x chinensis hybrids. Interestingly, despite repeated 
sampling, only a few hybrid individuals have been found in 
the native range (Gaskin, pers. comm.). Therefore, hybrid- 
ization appears to be a major part of the tamarisk invasion 
in the United States, but it is unknown exactly how hybrid- 
ization influences defense against herbivory (Gaskin and 
Kazmer 2009; Moran et al. 2009). In a set of experiments 
investigating cold hardiness and leaf phenology of hybrid 
tamarisk populations, Friedman et al. (2008, 2011) demon- 
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strated inherited latitudinal variation in these traits. When 
plants were grown in a common environment, southern 
populations had significantly increased susceptibility to 
cold temperatures and later onset of fall leaf senescence 
compared to northern populations. There is no indication 
to what degree this latitudinal variation results from species 
introgression or whether variation in these traits influences 
interactions with biological control (e.g., tolerance and re- 
growth following defoliation). 

We used a DNA fingerprinting technique to quantify 
species introgression among tamarisk plants collected from 
a wide geographical range in the western United States. 
Then, we compared introgression values against several 
environmental conditions at the collection sites for their 
power in predicting tamarisk performance traits when 
plants were grown in a common garden and subsequently 
subjected to defoliation experiments. Information on 
genetic differences in tamarisk growth and defense, along 
with knowledge of hybrid plant distribution, can poten- 
tially improve biological control in this system where evo- 
lution for at least some plant traits (i.e., cold hardiness, leaf 
phenology) has occurred. For instance, if tamarisk hybrid- 
ization occurs in a predictable pattern across the western 
United States and certain tamarisk hybrids have evolved 
more or less susceptibility to feeding by D. carinulata, this 
information can be applied to improve the efficiency of 
biological control tactics. 

Materials and methods 

Plant material 

A large common garden consisting of nearly 350 tamarisk 
shrubs located in Fort Collins, Colorado, was the source of 
plant material for this study (see Friedman et al. 2008, 
2011). All plants were originally collected from 15 natural 
populations along a latitudinal gradient from Montana to 
Texas where D. carinulata had yet to establish. Plants were 
3 years old at the onset of the current study. The large gar- 
den excluded vertebrate herbivores, and D. carinulata was 
never observed at the site. For the experiments described 
below, 72 genotypes representing 14 populations were ran- 
domly sampled from the large garden (Table SI). Six to ten 
cuttings (25 cm long, 2-8 mm diameter) were harvested 
from each of the selected genotypes. The cuttings were 
dipped in rooting hormone, planted individually in 100% 
perlite, and placed on a mist bench for 4 weeks to stimulate 
root growth. Survivors were transplanted into 13 
(W) x 13 (W) x 30 (H) cm pots with a mixture of 80:20 
of potting soil and sand. These plants were kept in the 
greenhouse under constant temperature and light for sev- 
eral weeks and then transferred to an outdoor shade house 
on the campus of Colorado State University, Fort Collins, 
Colorado. 
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AFLPs 

For DNA fingerprinting, we used clones of all available 
genotypes from the large common garden (n = 324; Fried- 
man et al. 2008) plus data available from other North 
American plants (n = 47) as well as native T. ramosissima 
(n = 56) and T. chinensis (n = 65) collected in Asia (Gas- 
kin and Kazmer 2009). Genomic DNA was extracted from 
approximately 20 mg of silica-dried plant material using a 
modified CTAB method (Hillis et al. 1996). Fragment PCR 
and AFLP analysis used the same protocol and platform 
found in Gaskin and Kazmer (2009). In short, primer pairs 
used for amplification were Msel + CTA/EcoRI + ACC 
and Msel + CTA/EcoRI + ACT, and fragments were ana- 
lyzed using an Applied Biosystems 3130 Genetic Analyzer 
(Foster City, CA, USA). GeneMapper v4.0 (Applied Biosys- 
tems) was used to visualize fragments for presence and 
length, while Structure v2.3.3 (Pritchard et al. 2000; Falush 
et al. 2003; Hubisz et al. 2009) was used to calculate assign- 
ment scores and introgression levels in terms of one of the 
two parent species, T. ramosissima (hereafter referred to as 
'introgression'). One can calculate T. chinensis introgres- 
sion by simple subtraction (100-%T. ramosissima intro- 
gression). Plants were considered hybrids if introgression 
levels were <0.9 or >0.1 and posterior probability intervals 
did not reach 1.0 (Pritchard et al. 2000; Blair and Hufbauer 
2010). 

Performance and tolerance experiment 

Ten populations representing a large geographical area 
were selected for the plant performance and tolerance 
experiment (Table SI). From these populations, 43 geno- 
types had at least three clones that survived the propagation 
technique outlined above. The three most vigorous clones 
per genotype were randomly assigned to one defoliation 
treatment: herbivore, chemical, and control. Finally, after 
clones were assigned treatments, they were further ran- 
domly assigned to experimental blocks housed in 100 
(W) x 100 (W) x 30 (H) cm wood box frames, which 
aided in treatment prescriptions and accounted for spatial 
heterogeneity in the outdoor shade house. There were nine 
total blocks (three per treatment). 

For the herbivore treatment, on July 16, 2010, we intro- 
duced 150 adult D. carinulata (number chosen based on 
our experience with the insect at field sites) into each of 
three blocks assigned to the herbivore defoliation treat- 
ment. The beetles were collected from an established field 
population near Palisade, Colorado, and were starved 24 h 
to minimize any effects of host preference based on prior 
feeding experience. To contain them within the blocks, we 
placed cages constructed from lightweight polyester cloth 
and plastic tubing over the box frames. We removed these 
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cages after 14 days when adults reached the end of their life 
cycle and larvae had not yet hatched from egg masses. Once 
emerged, larvae were free to move and feed among overlap- 
ping plant canopies within each block, while canopies from 
neighboring blocks were kept isolated. Larvae, which are 
not highly mobile, were never recorded outside of the 
blocks to which they were assigned. Blocks of the other 
defoliation treatments (chemical and control) were caged 
in the same manner and timing as those in the herbivore 
treatment. 

Because we anticipated that some clones in the herbivore 
treatment would receive little or no defoliation, we 
included a chemical treatment to measure plant response 
to nearly complete defoliation. For plants in the three 
blocks assigned to chemical defoliation treatment, the 
nonsystemic, defoliating herbicide, carfentrazone-ethyl 
(Aim™, FMC Corp. Philadelphia, PA, USA) was applied 
with a four-nozzle boom and a backpack sprayer using two 
passes at a volume of 187 L per hectare, which is consistent 
with field application of this herbicide (S. Nissen, pers. 
comm.). We temporarily moved the plants 200 m outside 
of the garden for approximately 2 h to avoid spray drift to 
nontarget plants. The timing of the application purpose- 
fully coincided with the height of beetle defoliation (August 
13, 2010). The remaining three blocks were assigned to the 
control treatment, and plants contained in these blocks did 
not receive any prescribed defoliation. 

From June to October 2010, measurements of damage 
and plant performance were recorded (see below). During 
this period, plants were watered every third day and sup- 
plemented with nutrients monthly using 20-20-20 NPK fer- 
tilizer. From October 2010 to April 2011, the wood boxes 
were filled with mulch to prevent roots from freezing, and 
water was supplemented as needed. In April 2011, the 
plants were uncovered, and the watering and fertilizing 
regime continued until the end of the experiment in June 
2011 when all plants were harvested. 

Quantifying damage, plant performance, and tolerance 

We used percentage of plant canopy damaged as a measure 
of defoliation. On September 3, 2010, we took two digital 
photographs (both in horizontal perspective) of each plant 
against a black background. The photographs were ana- 
lyzed using ImageJ, which quantifies pixilated area of digi- 
tal images based on color spectrum (vl.44, U.S. National 
Institute of Health, http://imagej.nih.gov/ij/). The whole 
canopy area was determined by setting the hue threshold 
from 21-102, while the green canopy area was measured 
with a hue threshold from 47-102. We averaged the whole 
and green canopy areas for pairs of pictures for each plant. 
By comparing the green canopy to the whole canopy area, 
we were able to assess the percentage of canopy damaged 
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for each plant, which included yellow and brown desiccated 
portions of plant canopy characteristic of D. carinulata 
damage. Our technique of quantifying defoliation matched 
closely to a subjective score of percent plant damage 
(r = 0.94, P < 0.0001). 

Annual canopy growth was used as a proxy for plant fit- 
ness, necessary for calculating tolerance to defoliation. We 
used canopy growth in lieu of reproductive output given 
the age of tree clones and the duration of this study. In 
June 2010 (before defoliation treatments were applied) and 
again in June 2011, plant height and the perpendicular can- 
opy widths were used to calculate canopy volume. This 
measure of canopy volume was strongly correlated with 
canopy area calculated using digital images and ImageJ 
software (r = 0.79, P < 0.0001). Canopy area was log- 
transformed in order to meet model assumptions of anova, 
and percent annual growth was calculated using these 
transformed scores of canopy size. 

We followed the suggestion of Strauss and Agrawal 
(1999) when determining plant tolerance to defoliation. 
Tolerance was quantified as the slope of the relationship 
between canopy damage and fitness measurement at 1 year 
postdefoliation treatment for each of the 43 genotypes. A 
separate slope was calculated for each genotype using three 
pairs of damage and associated fitness scores (i.e., herbi- 
vore, chemical, and control replicates). A positive slope 
indicates overcompensation for defoliation, while a nega- 
tive slope reveals undercompensation. 

To detect biomass allocation differences 1 year after 
defoliation treatments, we harvested all plants at the end of 
the experiment (June 2011). The roots were carefully 
washed of potting soil and placed in separate paper bags 
for drying. Green foliage and woody aboveground growth 
(excluding the original cutting) were also put in separate 
bags. The samples were placed in an oven at 55°C until they 
were dry (72 h). We used a digital balance to record the 
dry mass of green foliage, woody stems, coarse roots, and 
fine roots. Coarse roots were defined as those >1 mm 
diameter and fine roots as those <1 mm in diameter. Root- 
to-shoot ratios were calculated by dividing total below- 
ground biomass by total aboveground biomass. Cutting 
diameter only increased minimally over the course of the 
experiment and was not used to investigate the effects of 
defoliation or other factors. However, we included initial 
cutting diameter as covariate in models investigating over- 
all plant performance (see below). 

Resistance experiment 

We developed a second experiment to investigate resistance 
of tamarisk genotypes to herbivory by D. carinulata 
because of the difficulty in accurately measuring feeding 
damage by this particular herbivore. Four populations 
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from the large garden were chosen based upon their dis- 
persed geographical origin with 7-8 genotypes per popula- 
tion and up to eight clones per genotype surviving the 
propagation technique (Table SI). In August 2010, after 
these cuttings were transplanted and established in the out- 
door garden for 10 weeks, we clipped ~30 g of fresh green 
foliage from each plant, placed the material in a paper sack, 
and dried for 3 days at 40°C. Then, we powdered each 
sample using a coffee grinder, cleaning the instrument with 
70% EtOH and drying between samples. Next, we added 
5 g of dried plant material, 750 mg of agar, and 30 ml of 
H 2 0 to individual, flat-bottomed glass test tubes. After 
mixing the solution, the test tubes were placed in a steam 
bath for 20 min and then cooled to room temperature. 

Although our procedure likely destroyed or altered plant 
volatiles, it allowed us to produce uniform feeding media 
representing each individual clone grown in the garden. 
This media were used in an initial bioassay - the first of its 
kind for Diorhabda and Tamarix - to determine whether 
there were feeding differences among hybrid plants. Two 
identical pellets media were extracted from the tubes using 
a 5-mm punch. A single pellet was placed on top of filter 
paper inside of a 50-mm plastic Petri dish. One dish was 
randomly chosen as a control (no larvae), while the other 
received five, preweighed second-instar D. carinulata larvae 
collected from the field site in Palisade, Colorado. The lar- 
vae were starved for 24 h before the experiment began to 
minimize the effects of prior feeding experience. The Petri 
dishes were placed in an environmental chamber for 48 h 
with a daily temperature (29°C/21°C) and light cycle (14-h 
L/8-h D). After 48 h, larvae were removed and weighed. 
The pellets were then dried at 45°C for 3 day. The dry 
masses of the 'herbivore' and 'control' pellets were then 
compared to obtain percent consumed. Thus, we obtained 
both a resistance score and associated larval mass for sev- 
eral clones of each plant genotype. We developed this pellet 
bioassay in response to the difficulty in measuring actual 
damage to tamarisk plants by D. carinulata feeding (i.e., 
abscission of whole stems regardless of amount of feeding 
damage; difficulty of quantifying feeding damage on three- 
dimensional photosynthetic stems). Moreover, fresh- 
picked tamarisk shoots dry quickly (<1 h) and are not con- 
sumed by D. carinulata (pers. obs.); thus, our approach 
circumnavigated issues with using live plant material. 

Statistical analysis 

We used Akaike's information criterion (AIC) as a means 
of model selection to determine whether introgression suf- 
ficiently explained variation in plant traits compared to 
abiotic factors from the population origins (Burnham and 
Anderson 2002). We used PROC MIXED (SAS v9.3, SAS 
Institute, Cary, NC, USA) with maximum-likelihood esti- 
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mation of error to calculate AICc scores. We also calculated 
Akaike weights (wi) and evidence ratios, both of which are 
used as additional techniques for determining the most 
informative model (Burnham and Anderson 2002). Five 
model sets were considered for the model selection (Table 
S2). One model focused on species introgression, while the 
other four models investigated the effects of local adapta- 
tion to latitude (proxy of daylength), average number of 
annual frost-free days, average annual minimum tempera- 
ture, and elevation of sites where tamarisk was collected. 
To characterize temperature variation at the collection 
sites, we used the DAYMET interpolated daily low temper- 
ature data for the period 1980-2003 gridded at 1-km reso- 
lution (Thornton et al. 1997, http://www.daymet.org). 

After using model selection to determine whether 
hybridization or adaptation to abiotic conditions at site 
origins was more important in predicting plant growth 
traits, we then employed mixed-model analysis of covari- 
ance (ancova) for the investigation of defoliation damage 
and plant performance measurements using the statistical 
program JMP v9.0.2 (SAS Institute, Cary NC, USA). Defo- 
liation treatment and experimental block were considered 
fixed effects (i.e., both factors were tested across all possible 
levels), while plant subject was modeled as a random effect 
(i.e., levels sampled from a larger population). Diameter of 
stem cuttings at the onset of plant propagation was used as 
a covariate to control for the effect of initial size. In the 
resistance experiment, variables in models included the 
random effect of plant subject and either introgression or 
one of the four environmental variables. For tolerance, we 
used a weighted regression analysis to determine the 
strength of the relationship between species introgression 
and tolerance to defoliation. The weights were calculated as 
1/er 2 from each of the 43 tolerance regressions (damage vs. 
fitness for each genotype). Finally, we used Student's (--tests 
(a = 0.05) for post hoc comparisons among group means 
of treatment groups. 

Results 

AFLPs and environmental conditions at tamarisk sites 
Tamarisk introgression was highly correlated with latitude 
of plant origin (r = 0.87), annual frost-free days (r = 0.78), 
and minimum temperatures (r = 0.87; all P < 0.0001) and 
moderately correlated with elevation (r = 0.41; P < 0.01) 
at sites where plant material was collected. The environ- 
mental conditions were correlated with each other with the 
exception of elevation, which was only moderately 
correlated with annual frost-free days (r = 0.59, P = 0.03) 
and not correlated with latitude or minimum temperatures. 
The amount of tamarisk introgression increased 
with increasing latitude in the western United States (% 
T. ramosissima introgression = 4.6*latitude — 123; anova: 
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£1,340 = 987, P < 0.0001; Fig. 1, Table SI). Tamarisk from 
Colorado River, Texas (32.0°N), had the lowest amount of 
T. ramosissima introgression (22.7 ± 1.7%), while tama- 
risk from Musselshell River, Montana (46.4°N), had 
the highest amount of T. ramosissima introgression 
(92.9 ± 1.1%). Mean introgression across all latitudes was 
59.3 ± 1.4%. Out of the 342 U.S. tamarisk shrubs included 
in the AFLP analysis, there were 37 individuals with 
assignment scores >0.9 for T. ramosissima (all from sites 
41.3-47.6°N) and only four individuals with assignment 
scores >0.9 for T. chinensis (all from 32.0-35.5 °N). 

Model selection 

Because introgression was highly correlated with latitude, 
frost-free days, and minimum temperatures of population 
origin, these predictor variables could not be placed in the 
same model investigating plant performance, thus justify- 
ing the use model selection. AIC scores showed that 
Model 1 (introgression) was more informative than the 
other models for root-to-shoot ratio, resistance, and larval 
mass (Table S3). For plant biomass, introgression and lati- 
tude were essentially equal in their ability to fit the 
observed data (A; AIC < 2). Only in one case (canopy 
growth rate) was an environmental variable, elevation, 
more informative than introgression. All models predict- 
ing defoliation damage had A, AIC < 2, indicating no 
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Figure 1 The relationship between tamarisk species introgression, 
expressed as percentage T. ramosissima, arid latitude of population ori- 
gin (%T. ramosissima introgression = 4.6*latitude - 123; R 2 = 0.74; 
n = 342). See Table S1 for detailed site descriptions. 



clear distinction among these candidates (Burnham and 
Anderson 2002). In summary, introgression explained as 
much or more variation in plant response to defoliation 
than abiotic factors from population origins; thus, we 
present the reminder of the results in terms of tamarisk 
introgression. 

Plant performance traits 

We provide full ancova reports for all measured plant 
responses in Appendix SI. Initial plant size (cutting diame- 
ter) had a significant influence on all aspects of plant bio- 
mass, while it had no effect on defoliation damage, growth 
rate, or root-to-shoot ratio. The random effect of plant 
subject explained 20-60% of variation across all plant 
traits. 

The introgression*treatment interaction was not signifi- 
cant for defoliation damage, indicating no difference in 
damage across introgression levels, nor was there a signifi- 
cant effect of introgression alone on this response variable 
(Appendix SI). However, there was a significant effect of 
treatment type on plant damage (£2,77.6 = 120, 
P < 0.0001). Damage from the chemical treatment (LS 
mean ± SE: 63.9 ± 2.6%, range: 20-93%) was greater 
than damage from the herbivore treatment (LS 
mean ± SE: 35.1 ± 2.6%, range: 3-94%). Both chemical 
and herbivore treatments resulted in significantly greater 
plant damage than the control treatment (LS mean ± SE: 
7.5 ± 2.6%, range: 1-32%). 

Our proxy for plant fitness, annual canopy growth rate, 
was also not affected by species introgression nor was there 
a significant interaction between introgression and defolia- 
tion (Appendix SI). Defoliation treatment had a significant 
effect on growth rate (£2,77.6 = 120, P < 0.0001). Growth 
rates between the treatments were 37.7 ± 5.0%, 
33.3 ± 5.0%, and 11.63 ± 6.1%, for control, herbivore, 
and chemical defoliation, respectively. Fitness scores ranged 
between —100.0 and 93.6% annual canopy growth. Six of 
the 43 clones assigned to the chemical treatment died dur- 
ing the winter following the treatment application (fit- 
ness = —100% canopy growth). These plants were retained 
in the analysis of tolerance to defoliation (see below). 

Total plant biomass decreased with species introgres- 
sion (biomass = 21.0-0. 17*introgression; £1,39.7 = 4.7, 
P = 0.04). In addition, defoliation treatment had a signifi- 
cant effect on plant biomass (F 2 , 74.5 = 16.6, P < 0.0001), 
but the interaction between treatment and introgression 
was not significant. Defoliation by either chemical or herbi- 
vore treatments led to significantly reduced coarse root and 
woody stem biomass (Fig. 2). On the other hand, only 
plants that were chemically defoliated had significantly 
lower green foliage and fine root biomass compared to 
control plants. 
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Green foliage Woody stems Coarse roots Fine roots 

Figure 2 LS mean biomass ± SE of four plant tissue groups for three 
defoliation treatments. Chemical: black bars; herbivore: dark gray bars; 
control: light gray bars. Different letters correspond to differences 
between means within a plant tissue category (Student's t-test, 
a. = 0.05). Defoliation treatments were applied in the summer of 2010, 
and plant parts were harvested, dried, and weighed the following 
spring. Mean percent canopy defoliation was 64%, 35%, and 7.5% for 
chemical, herbivore, and control treatments, respectively. 

Results of ancova exposed both a significant treatment 
effect (F 2 ,78.6 = 8.5, P < 0.001) and a significant effect of 
species introgression (F lA03 = 4.5, P = 0.04) on root-to- 
shoot ratio, but not a significant interaction between these 
two factors. Plants in the herbivore (LS mean: 0.29 ± 0.02) 
and control (LS mean: 0.27 ± 0.02) treatment groups both 
had significantly greater root-to-shoot ratios than plants in 
the chemical treatment (LS mean: 0.19 ± 0.02), but were 
not significantly different than each other. Introgression 
showed a positive relationship with root-to-shoot ratios 
(ratio = 0.001 16*introgression + 0.20), implying that 
plants with high levels of T. ramosissima introgression 
invest more in belowground rather than in aboveground 
growth. 

Tolerance 

We provide a summary of tolerance slopes and standard 
errors for each of the 43 genotypes in Appendix S2. These 
genotypes, which were haphazardly chosen from the larger 
pool of 342 plants, had a mean species introgression level 
of 58.0 ± 3.4% T. ramosissima (range: 20.5-96.2%). Toler- 
ance demonstrated a positive relationship with introgres- 
sion using weighted regression analysis (tolerance = 0.004* 
introgression - 0.40; F 1>41 = 6.27, P = 0.02). Plants with 



high levels of T. ramosissima introgression had increased 
tolerance to defoliation (Fig. 3). 

Resistance 

High resistance to herbivory was reflected by a small pro- 
portion, p, of plant material consumed in the pellet bioas- 
say. There was a negative relationship between species 
introgression and resistance (fi,28.9 = 15.2, P < 0.001; 
Fig. 4A): Pellets with high levels of T. ramosissima were 
more completely consumed. Additionally, larvae used in 
this experiment gained more mass when fed pellets with 
high T. ramosissima introgression (-Fi,28.9 = 8.0, P < 0.01; 
Fig. 4B), confirming the low resistance of these plant geno- 
types. 

Discussion 

Tamarisk evolution via hybridization is widespread in 
North America. Using markers for a single-locus nuclear 
DNA gene (PepC), Gaskin and Schaal (2002) first reported 
that the most common plant in the U.S. invasion was a 
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Figure 3 Tolerance to defoliation across species introgression for plant 
genotypes in the experiment. Positive tolerance values correspond to 
plant overcompensation, while negative values correspond to plant un- 
dercompensation to defoliation. Tolerance demonstrated a positive 
relationship with introgression when the SEs of tolerance scores were 
included in a weighted regression (tolerance = 0.004*introgression - 
0.40; R 2 = 0.13). Weighted regression was used because some clones 
in the chemical treatment died, thus leading to more error for those 
particular genotypes (triangles). See Appendix S2 for detailed tolerance 
scores of each genotype. 
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Figure 4 Plant genotype resistance to herbivory by D. carinulata in the 
pellet bioassay (A). Resistance is defined as 1 -proportion of pellet 
damage. Pellets made from plants with high levels of T. ramosissima 
introgression had greater resistance (resistance = 0.91 9-0. 002*intro- 
gression; F, 2 s.9 = 1 5.2, P < 0.001; R 2 = 0.42). Some resistance values 
were >1 in some cases because the treatment pellet weighed slightly 
more than the control pellet and it was not consumed. Percent 
change in larvae mass in the pellet bioassay (B). Plant species introgres- 
sion had a negative effect on larval performance (%change in 
larval mass = 0.1 5*introgression - 15.2; F 1|27 . 8 = 8.0, P<0.01; 
R 2 = 0.29). 

hybrid of two species-specific genotypes that were geo- 
graphically isolated in the native range of Asia. Then, as fol- 
low-up, Gaskin and Kazmer (2009) used a multilocus 



approach across all genomic DNA (i.e., AFLP) and Bayes- 
ian cluster analysis to determine that 83-87% of the U.S. 
tamarisk genotypes are indeed hybrids between X. ramo- 
sissima and X. chinensis parent species. In the present 
study, which used the same AFLP protocol and laboratory 
as Gaskin and Kazmer (2009), we found that 89% of the 
371 North American plants in the study were hybrids with 
only handful of individuals (10% and 1%, respectively) 
classified as X. ramosissima and X. chinensis parent species. 
Additionally, we found a significant relationship between 
latitude of tamarisk populations and interspecific hybrid- 
ization, following similar studies (Friedman et al. 2008; 
Gaskin and Kazmer 2009); however, we describe the pat- 
tern in detail for the first time: Among tamarisk popula- 
tions in the western United States, every one degree 
increase in latitude between 32.0 and 46.7°N corresponds 
to a 4.6% positive change in X. ramosissima introgression 
(Fig. 1). As our experiments revealed, the relationship 
between latitude and tamarisk hybridization has direct 
impact on interactions with a specialized herbivore intro- 
duced for control of this invasive plant. 

The exact cause of the latitudinal cline in tamarisk 
hybridization is unknown. One possible mechanism is that 
hybridization occurred in nurseries in the early 1800s 
before tamarisk was introduced at multiple locations in the 
western United States (Gaskin et al. 2012), while another is 
that several independent hybridization events could have 
occurred in areas where both parent species were co-intro- 
duced (Gaskin and Kazmer 2009; Friedman et al. 2011). 
Under either scenario, subsequent sorting out of adapted 
hybrid genotypes then occurred postintroduction. The 
abundance of tamarisk hybrids in North America com- 
pared to the lack of interspecific hybridization in the native 
range in Asia suggests that species introgression has played 
an important role for invasion success (Gaskin and Kazmer 
2009; Gaskin et al. 2012). Similar patterns of latitude and 
genetic variation are demonstrated by at least one other 
invasive shrub, Eupatorium adenophorum (Huang et al. 
2009). 

We used Akaike's information criterion (AIC) in a 
model selection approach to determine the relative 
strengths of the correlated factors of introgression and abi- 
otic conditions at population origins (i.e., hybridization 
and local adaptation) in explaining plant growth and 
defense traits (Burnham and Anderson 2002). Introgression 
was more informative than other models in many cases, 
including those for tissue allocation and resistance (Table 
S3). On the other hand, growth rate was strongly related to 
elevation, which most likely resulted from the unimodal 
relationship between elevation and latitude of collection 
sites. For total biomass, introgression and latitude 
explained equal amounts of variation and both had a sig- 
nificant negative relationship with this plant trait (Table 
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S3). Latitudinal clines in similar traits have been recorded 
in a number of invasive plant species in their introduced 
ranges (Weber and Schmid 1998; Kollmann and Baiiuelos 
2004; Maron et al. 2004; Leger and Rice 2007; Montague 
et al. 2008; Keller et al. 2009). We propose that in the case 
of tamarisk latitudinal clines, interspecific hybridization 
may have provided the sufficient genetic material for selec- 
tion to act upon following introduction to North America. 
Regardless of the tight correlation between latitude and 
introgression, genetic information (e.g., AFLPs), and not 
simply geographical location, was more informative for 
root-to-shoot ratios and herbivore defense strategies for 
hybrid tamarisk plants. 

In the outdoor garden experiment, there were a signifi- 
cant decrease in total plant biomass and a significant 
increase in root-to-shoot ratios in plants grown with 
increasing percentage of T. ramosissima alleles. Sexton 
et al. (2002) observed the variation in root-to-shoot ratios 
of comparable tamarisk populations, with northern 
populations investing more in belowground growth. Addi- 
tionally, we found that plants with high levels of T. ramo- 
sissima introgression were more tolerant of defoliation 
than plants with low levels of T. ramosissima introgression 
(Fig. 3). Given that belowground carbohydrate stores are 
important for tamarisk leaf flush following defoliation 
(Hudgeons et al. 2007), tamarisk hybridization is impor- 
tant in an ecological and evolutionary sense. For instance, 
plants from northern populations may experience more 
dieback as the result of extreme low temperatures, 
selecting for genotypes with greater belowground alloca- 
tion and root-to-shoot ratios (Friedman et al. 2008). 
Indeed, T. ramosissima introgression was associated with 
higher latitudes and higher root-to-shoot ratios in our 
experiment. T. ramosissima genotypes may have adapted 
strategies to cope with D. carinulata outbreaks in the 
native range where they co-occur (Lewis et al. 2003), 
while in the introduced range, D. carinulata has caused 
the most impact on tamarisk populations whose introgres- 
sion values are <60% T. ramosissima (DeLoach et al. 
2009; Gaskin and Kazmer 2009). We have begun investi- 
gating levels of introgression at field sites where D. carinu- 
lata has failed to establish to determine whether this 
pattern holds true. In addition, at least five species of Dio- 
rhabda co-occur with numerous Tamarix spp. in the 
native range of Asia and Europe, but no Diorhabda species 
has been recorded on T. chinensis in eastern Asia (Tracy 
and Robbins 2009). Thus, T. c/imenszs-dominated 
genotypes may invest more in chemical defenses against 
herbivores in general than to allocating resources to 
regrowth following aboveground dieback. Given that tol- 
erance is an herbivore defense strategy that has a strong 
degree of inheritance (Strauss and Agrawal 1999; Agrawal 
et al. 2004), hybridization between T. ramosissima and 
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T. chinensis could affect the outcome of biological control 
at certain field sites. 

In the resistance experiment, the overall trend in larval 
mass indicated that, in fact, most insects lost weight during 
the 48-h experiment, indicating that plant pellets would be 
a poor choice for an artificial diet. Still, differences in 
weight gained (or lost) were evident among plant geno- 
types, and these differences coincided with tamarisk intro- 
gression. Plant pellets with high levels of T. ramosissima 
introgression were consumed in higher quantities than their 
counterparts (Fig. 4A). Similarly, larvae that were fed pel- 
lets made from plants with high T. ramosissima introgres- 
sion performed better (Fig. 4B). Several accessions and 
species of tamarisk have been tested in D. carinulata feeding 
trials, oftentimes with significant differences in attack or 
insect development (DeLoach et al. 2003; Lewis et al. 2003; 
Milbrath and DeLoach 2006; Dalin et al. 2009; Moran et al. 
2009). We did not detect differences in resistance among 
genotypes for the outdoor study, even though we used a 
precise digital imaging tool to objectively assess damage. 
One explanation is that when tissue is damaged by D. cari- 
nulata feeding, a common plant response is to abscise entire 
leaves and stems. Thus, beetles in the outdoor experiment 
could have preferentially fed upon various tamarisk 
hybrids, but the amount of herbivory was masked by this 
abscission. Nonetheless, the bioassay experiment confirmed 
differences in resistance among hybrid genotypes, and these 
differences were related to tamarisk introgression. 

Differences in herbivore resistance among hybrid plant 
genotypes occur in natural systems (Fritz et al. 1998; 
McGuire and Johnson 2006; Krebs et al. 2011). A possible 
mechanism for these differences is novel chemistry of 
hybrid host plants (Orians 2000). Three reviews involving 
over 80 plant taxa show that 5.5% of secondary compounds 
of hybrid genotypes are novel compared to either parent 
species, while nearly 20% of hybrid plants either over- or 
underexpress secondary compounds compared to their 
parents. Moreover, roughly 5-6% of hybrid genotypes have 
higher resistance to herbivory compared to parent species 
(Orians 2000; Cheng et al. 2011). Because herbivorous 
insects use plant chemistry as a means of host recognition, 
novel chemical combinations due to hybridization may 
make plants unrecognizable or unpalatable to highly spe- 
cialized herbivores (Schoonhoven et al. 2005). While plant- 
based chemical attractants are known for D. carinulata 
(Cosse et al. 2006), there is no information whether these 
chemicals differ among hybrid genotypes. Furthermore, 
there are no published records of tamarisk chemistry rela- 
tive to insect herbivory in the introduced range. It would 
be interesting to see whether differences in tamarisk resis- 
tance to herbivory can be explained by differences in the 
production of important secondary metabolites by hybrid 
tamarisk plants. 
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Finally, our study also revealed an interesting ecological 
relationship between defoliation and plant performance, 
which may have implications in an applied sense. Only a 
modest amount of damage by beetles (mean = 35.1%) sig- 
nificantly reduced annual growth of woody stems and 
coarse roots compared to control plants (mean dam- 
age = 7.5%, most likely attributed to grasshoppers observed 
in the outdoor shade house). Moreover, beetle defoliation 
had the same effect as did a larger amount of damage 
(mean = 63.9%) by chemical defoliation (Fig. 2). There- 
fore, D. carinulata appears to reduce allocation to these 
structural and storage tissues as much as chemical 
defoliation. On the other hand, chemical defoliation by 
carfentrazone-ethyl significantly reduced the growth of 
photosynthetic tissue and fine roots the following year, 
while defoliation by herbivores did not. While chemical 
control appears to have a more immediate effect in terms of 
overall plant size, our experiment did not take into account 
that plants often experience repeated herbivore defoliation 
in areas where the insect is established. At field sites where 
D. carinulata has had the most impact, the insect completes 
2-4 generations per year while depleting carbohydrate 
reserves through repeated defoliation, ultimately leading to 
plant mortality in as little as 3 years (Lewis et al. 2003; 
Hudgeons et al. 2007; Pattison et al. 2011). 

Adaptation appears to be an important evolutionary 
mechanism for hybrid tamarisk populations in terms of 
particular life history traits (e.g., cold hardiness and leaf 
phenology; Friedman et al. 2008, 2011) and tolerance and 
resistance to a specialized herbivore (this study). The 
source of genetic variation for selection to act upon in this 
system may have been due to interspecific hybridization. 
Quantitative trait loci mapping and breeding experiments 
involving parent species to produce Fl, F2, and back- 
crossed offspring are necessary to be definitive (Cheng 
et al. 2011). Additionally, comparing parent species from 
the native range to parental and hybrid genotypes from the 
introduced range would reveal how important hybridiza- 
tion is in this system (Keller and Taylor 2008; Ridley and 
Ellstrand 2010). In terms of assessing the strength of latitu- 
dinal clines, caution should be taken when drawing infer- 
ences from a single common garden experiment (Williams 
et al. 2008; Colautti et al. 2009). Multiple common gardens 
with plants collected from a wide range of latitudes in both 
ranges are ideal for determining whether evolution of lati- 
tudinal clines has occurred (Maron et al. 2004; Colautti 
et al. 2009). Given that a study which combines both lati- 
tude and hybridization may be difficult, perhaps focusing 
attention on U.S. populations that naturally possess a wide 
range of genetic variability could be fruitful. 

Rapid evolution can affect outcomes in the application 
of biological control (see July 2012 special issue of Evolu- 
tionary Applications). Curiously, the tamarisk system has 



been an example of two such cases, both involving the biol- 
ogy of the insect (Thomas et al. 2010; Bean et al. 2012). 
Here, we show that evolution of tamarisk via hybridization 
affects interactions between plants and a specialized herbi- 
vore introduced for its control. High tolerance of northern 
tamarisk (i.e., T. ramosissima-dommated genotypes) is 
apparently related to belowground tissue allocation, while 
high herbivore resistance of southern tamarisk (i.e., 
T. chinensis-dominated genotypes) may be due to chemical 
or structural defenses. Intense herbivory by D. carinulata 
may cause rapid natural selection of hybrid plants combin- 
ing tolerance and resistance more successfully than either 
parent species alone. In the short term, if hybrids are 
spread across the landscape in a predictable manner, as is 
the case with Tamarix invasion in the western United 
States, this knowledge regarding the response of particular 
hybrids to herbivory can be exploited to maximize biologi- 
cal control efforts. 

Acknowledgements 

Funding was provided by the Colorado State University 
Agricultural Experiment Station and the NSF Research 
Coordination Network. We are grateful to R. Hufbauer, P. 
McEvoy, E. Coombs, F. Grevstad, and three anonymous 
reviewers who provided comments on the manuscript. D. 
Bean, L. Bjostad, and S. Nissen assisted with supplying 
insects, providing bioassay materials, and applying herbi- 
cide, respectively. 

Data archiving statement 

Data for this study are available from the Dryad Digital 
Repository: http://dx.doi.org/10.5061/dryad.vc62h. 

Literature cited 

Abbott, R. J., J. K. James, R. I. Milne, and A. C. M. Gillies 2003. Plant 
introductions, hybridization and gene flow. Philosophical Transac- 
tions of the Royal Society of London Series B-Biological Sciences 
358:1123-1132. 

Agrawal, A. A., J. K. Conner, and J. R. Stinchcombe 2004. Evolution of 
plant resistance and tolerance to frost damage. Ecology Letters 
7:1199-1208. 

Bean, D. W., P. Dalin, and T. L. Dudley 2012. Evolution of critical day 
length for diapause induction enables range expansion of Diorhabda 
carinulata, a biological control agent against tamarisk (Tamarix spp.). 
Evolutionary Applications 5:51 1—523. 

Blair, A. C, and R. A. Hufbauer 2010. Hybridization and invasion: one 
of North America's most devastating invasive plants shows evidence 
for a history of interspecific hybridization. Evolutionary Applications 
3:40-51. 

Blair, A. C, U. Schaffner, P. Hafliger, S. K. Meyer, and R. A. Hufbauer 
2008. How do biological control and hybridization affect enemy 
escape? Biological Control 46:358-370. 



390 



©201 4 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 381-393 



Williams et al. 

Blossey, B., and R. Notzold 1995. Evolution of increased competitive 
ability in invasive nonindigenous plants - a hypothesis. Journal of 
Ecology 83:887-889. 

Bossdorf, O., H. Auge, L. Lafuma, W. E. Rogers, E. Siemann, and D. 
Prati 2005. Phenotypic and genetic differentiation between native and 
introduced plant populations. Oecologia 144:1—11. 

Burnham, K. P., and D. R. Anderson 2002. Model Selection and Multi- 
model Inference. 2nd edn. Springer- Verlag, New York. 

Campbell, L. G., A. A. Snow, and C. E. Ridley 2006. Weed evolution after 
crop gene introgression: greater survival and fecundity of hybrids in a 
new environment. Ecology Letters 9:1198—1209. 

Cheng, D. D., K. Vrieling, and P. G. L. Klinkhamer 2011. The effect of 
hybridization on secondary metabolites and herbivore resistance: 
implications for the evolution of chemical diversity in plants. Phyto- 
chemistry Reviews 10:107—117. 

Clausen, J. D., D. Keck, and W. M. Hiesey 1940. Experimental Studies on 
the Nature of Species I. The Effects of Varied Environments on North 
American Plants. Carnegie Institute of Washington, Washington, DC. 

Colautti, R. I., J. L. Maron, and S. C. H. Barrett 2009. Common garden 
comparisons of native and introduced plant populations: latitudinal 
clines can obscure evolutionary inferences. Evolutionary Applications 
2:187-199. 

Cosse, A. A., R. J. Bartelt, B. W. Zilkowski, D. W. Bean, and E. R. An- 
dress 2006. Behaviorally active green leaf volatiles for monitoring the 
leaf beetle, Diorhabda elongata, a biocontrol agent of saltcedar, Ta- 
marix spp. Journal of Chemical Ecology 32:2695—2708. 

Cuda, J. P., L. R. Christ, V. Manrique, W. A. Overholt, G. Wheeler, and 
D. A. Williams 2012. Role of molecular genetics in identifying 'fine 
tuned' natural enemies of the invasive Brazilian peppertree, Schinus 
terebinthifolius: a review. BioControl 57:227—233. 

Dalin, P., M. J. O'Neal, T. Dudley, and D. W. Bean 2009. Host plant 
quality of Tamarix ramosissima and T. parviflora for three sibling spe- 
cies of the biocontrol insect Diorhabda elongata (Coleoptera: Chryso- 
melidae). Environmental Entomology 38:1373—1378. 

DeBach, P. E. 1964. Biological Control of Insect Pests and Weeds. Rein- 
hold Publishing Corporation, New York. 

DeLoach, C. J., P. A. Lewis, J. C. Herr, R. I. Carruthers, J. L. Tracy, and J. 
Johnson 2003. Host specificity of the leaf beetle, Diorhabda elongata 
deserticola (Coleoptera: Chrysomelidae) from Asia, a biological con- 
trol agent for saltcedars {Tamarix: Tamaricaceae) in the Western Uni- 
ted States. Biological Control 27:117-147. 

DeLoach, C. J., A. E. Knutson, and P. J. Moran. 2009. Progress on Bio- 
logical Control of Saltcedar. USDA and Texas AgriLife, http://www. 
orangentia.net/SaltcedarPublicReportHandoutll-10-09.pdf (accessed 
on July 30 2013). 

Donovan, L. A., D. M. Rosenthal, M. Sanchez- Velenosi, L. H. Rieseberg, 
and F. Ludwig 2010. Are hybrid species more fit than ancestral parent 
species in the current hybrid species habitats? Journal of Evolutionary 
Biology 23:805-816. 

Dudley, T. L., D. W. Bean, R. R. Pattison, and A. Caires 2012. Selectivity 
of a biological control agent, Diorhabda carinulata Desrochers, 1870 
(Coleoptera: Chrysomelidae) for host species within the genus Ta- 
marix Linneaus, 1753. Pan-Pacific Entomologist 88:319-341. 

Ellstrand, N. C, and K. A. Schierenbeck 2000. Hybridization as a stimu- 
lus for the evolution of invasiveness in plants? Proceedings of the 
National Academy of Sciences of the United States of America 
97:7043-7050. 

Falush, D., M. Stephens, and J. K. Pritchard 2003. Inference of popula- 
tion structure using multilocus genotype data: linked loci and corre- 
lated allele frequencies. Genetics 164:1567—1587. 



Invasive plant hybridization and biological control 

Friedman, J. M., G. T. Auble, P. B. Shafroth, M. L. Scott, M. F. Merigli- 
ano, M. D. Preehling, and E. K. Griffin 2005. Dominance of non- 
native riparian trees in western USA. Biological Invasions 7:747—751. 

Friedman, J. M., J. E. Roelle, J. F. Gaskin, A. E. Pepper, and J. R. Manhart 

2008. Latitudinal variation in cold hardiness in introduced Tamarix 
and native Populus. Evolutionary Applications 1:598—607. 

Friedman, J. M., J. E. Roelle, and B. S. Cade 2011. Genetic and environ- 
mental influences on leaf phenology and cold hardiness of native and 
introduced riparian trees. International Journal of Biometeorology 
55:775-787. 

Fritz, R. S., B. M. Roche, and S. J. Brunsfeld 1998. Genetic variation in 
resistance of hybrid willows to herbivores. Oikos 83:117-128. 

Fritz, R. S., C. Moulia, and G. Newcombe 1999. Resistance of hybrid 
plants and animals to herbivores, pathogens, and parasites. Annual 
Review of Ecology and Systematics 30:565-591. 

Gaskin, J. F., and D. J. Kazmer 2009. Introgression between invasive salt- 
cedars {Tamarix chinensis and T. ramosissima) in the USA. Biological 
Invasions 11:1121-1130. 

Gaskin, J. F., and B. A. Schaal 2002. Hybrid Tamarix widespread in US 
invasion and undetected in native Asian range. Proceedings of the 
National Academy of Sciences of the United States of America 
99:11256-11259. 

Gaskin, J. F., A. S. Birken, and D. J. Cooper 2012. Levels of novel hybrid- 
ization in the saltcedar invasion compared over seven decades. Biolog- 
ical Invasions 14:693-699. 

Hillis, D. M., B. K. Mable, A. Larson, S. K. Davis, and E. A. Zimmer. 
1996. Nucleic acids IV: sequencing and cloning. In D. M. Hillis, C. 
Moritz, and B. K. Mable, eds. Molecular Systematics, pp. 321-381. 
Sinauer Associates, Sunderland, MA. 

Hodgins, K. A., and L. Rieseberg 2011. Genetic differentiation in life-his- 
tory traits of introduced and native common ragweed {Ambrosia ar- 
temisiifolia) populations. Journal of Evolutionary Biology 24:2731- 
2749. 

Huang, W. K., F. H. Wan, J. Y. Guo, B. D. Gao, B. Y. Xie, and D. L. Peng 

2009. AFLP analyses of genetic variation of Eupatorium adenophorum 
(Asteraceae) populations in China. Canadian Journal of Plant Science 
89:119-126. 

Hubisz, M. J., D. Falush, M. Stephens, and J. K. Pritchard 2009. Inferring 
weak population structure with the assistance of sample group infor- 
mation. Molecular Ecology Resources 9:1322—1332. 

Hudgeons, J. L., A. E. Knutson, K. M. Heinz, C. J. DeLoach, T. L. Dud- 
ley, R. R. Pattison, and J. R. Kiniry 2007. Defoliation by introduced 
Diorhabda elongata leaf beetles (Coleoptera: Chrysomelidae) reduces 
carbohydrate reserves and regrowth of Tamarix (Tamaricaceae). Bio- 
logical Control 43:213-221. 

Hull-Sanders, H. M., R. Clare, R. H. Johnson, and G. A. Meyer 2007. 
Evaluation of the evolution of increased competitive ability (EICA) 
hypothesis: loss of defense against generalist but not specialist herbi- 
vores. Journal of Chemical Ecology 33:781—799. 

Hultine, K. R., T. L. Dudley, and S. W. Leavitt 2013. Herbivory- induced 
mortality increases with radial growth in an invasive riparian phreato- 
phyte. Annals of Botany 111:1197-1206. 

Keller, S. R., and D. R. Taylor 2008. History, chance and adaptation dur- 
ing biological invasion: separating stochastic phenotypic evolution 
from response to selection. Ecology Letters 11:852—866. 

Keller, S. R., D. R. Sowell, M. Neiman, L. M. Wolfe, and D. R. Taylor 
2009. Adaptation and colonization history affect the evolution of 
clines in two introduced species. New Phytologist 183:678-690. 

van Klinken, R. D., and S. Raghu 2006. A scientific approach to agent 
selection. Australian Journal of Entomology 45:253-258. 



©201 4 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 381-393 



391 



Invasive plant hybridization and biological control 



Williams et al. 



Kollmann, J., and M. J. Banuelos 2004. Latitudinal trends in growth and 

phenology of the invasive alien plant Impatiens glandulifera (Balsam- 

inaceae). Diversity and Distributions 10:377—385. 
Krebs, C„ E. Gerber, D. Matthies, and U. Schaffner 2011. Herbivore 

resistance of invasive Fallopia species and their hybrids. Oecologia 

167:1041-1052. 

Lee, C. E. 2002. Evolutionary genetics of invasive species. Trends in Ecol- 
ogy 8c Evolution 17:386-391. 

Leger, E. A., and K. J. Rice 2007. Assessing the speed and predictability 
of local adaptation in invasive California poppies {Eschscholzia califor- 
nica). Journal of Evolutionary Biology 20:1090—1103. 

Lewis, P. A., C. J. DeLoach, A. E. Knutson, J. L. Tracy, and T. O. Robbins 

2003. Biology of Diorhabda elongate deserticola (Coleoptera: Chryso- 
melidae), an Asian leaf beetle for biological control of saltcedars 
(Tamarix spp.) in the United States. Biological Control 27:101-116. 

Mack, R. N., D. Simberloff, W. M. Lonsdale, H. Evans, M. Clout, and F. 
A. Bazzaz 2000. Biotic invasions: causes, epidemiology, global conse- 
quences, and control. Ecological Applications 10:689-710. 

Maron, J. L., M. Vila, R. Bommarco, S. Elmendorf, and P. Beardsley 

2004. Rapid evolution of an invasive plant. Ecological Monographs 
74:261-280. 

McGuire, R. J., and M. T. J. Johnson 2006. Plant genotype and induced 
responses affect resistance to herbivores on evening primrose {Oeno- 
thera biennis). Ecological Entomology 31:20—31. 

Milbrath, L. R., and C. J. DeLoach 2006. Host specificity of different 
populations of the leaf beetle Diorhabda elongata (Coleoptera: Chryso- 
melidae), a biological control agent of saltcedar {Tamarix spp.). Bio- 
logical Control 36:32-48. 

Montague, J. L., S. C. H. Barrett, and C. G. Eckert 2008. Re-establish- 
ment of clinal variation in flowering time among introduced popula- 
tions of purple loosestrife {Lythrum salicaria, Lythraceae). Journal of 
Evolutionary Biology 21:234-245. 

Moran, P. J., C. J. DeLoach, T. L. Dudley, and J. Sanabria 2009. Open 
field host selection and behavior by tamarisk beetles {Diorhabda spp.) 
(Coleoptera: Chrysomelidae) in biological control of exotic saltcedars 
{Tamarix spp.) and risks to non-target athel {T. aphylla) and native 
Frankenia spp. Biological Control 50:243—261. 

Muller-Scharer, H., U. Schaffner, and T. Steinger 2004. Evolution in 
invasive plants: implications for biological control. Trends in Ecology 
& Evolution 19:417^22. 

Nagler, P. L., E. P. Glenn, C. S. Jarnevich, and P. B. Shafroth. 2011. Dis- 
tribution and abundance of saltcedar and Russian olive in the western 
United States. Critical Reviews in Plant Sciences 30:508—523. 

Orians, C. M. 2000. The effects of hybridization in plants on secondary 
chemistiy: Implications for the ecology and evolution of plant-herbi- 
vore interactions. American Journal of Botany 87:1749—1756. 

Pattison, R. R., C. M. D'Antonio, and T. L. Dudley 2011. Biological con- 
trol reduces growth, and alters water relations of the saltcedar tree 
{Tamarix spp.) in western Nevada, USA. Journal of Arid Environ- 
ments 75:346-352. 

Pimentel, D., R. Zuniga, and D. Morrison 2005. Update on the environ- 
mental and economic costs associated with alien -invasive species in 
the United States. Ecological Economics 52:273-288. 

Prentis, P. J., J. R. U. Wilson, E. E. Dormontt, D. M. Richardson, and A. 
J. Lowe 2008. Adaptive evolution in invasive species. Trends in Plant 
Science 13:288-294. 

Pritchard, J. K., M. Stephens, and P. Donnelly 2000. Inference of popula- 
tion structure using multilocus genotype data. Genetics 155:945-959. 

Ridley, C. E., and N. C. Ellstrand 2010. Rapid evolution of morphology 
and adaptive life history in the invasive California wild radish {Raph- 



anus sativus) and the implications for management. Evolutionary 

Applications 3:64-76. 
Rieseberg, L. H., S. C. Kim, R. A. Randell, K. D. Whitney, B. L. Gross, C. 

Lexer, and K. Clay 2007. Hybridization and the colonization of novel 

habitats by annual sunflowers. Genetica 129:149—165. 
Robinson, T. W. 1965. Introduction, spread and areal extent of saltcedar 

{Tamarix) in the western states. US Geological Survey Professional 

Paper 491 -A. 

Schierenbeck, K. A., and N. C. Ellstrand 2009. Hybridization and the 
evolution of invasiveness in plants and other organisms. Biological 
Invasions 11:1093-1105. 

Schoonhoven, L. M., J. J. A. van Loon, and M. Dicke 2005. Insect-Plant 
Biology. Oxford University Press, Oxford. 

Sexton, J. P., J. K. McKay, and A. Sala 2002. Plasticity and genetic diver- 
sity may allow saltcedar to invade cold climates in North America. 
Ecological Applications 12:1652-1660. 

Snow, A. A., B. Andersen, and R. B. Jorgensen 1999. Costs of transgenic 
herbicide resistance introgressed from Brassica napus into weedy B. 
rapa. Molecular Ecology 8:605-615. 

Strauss, S. Y., and A. A. Agrawal 1999. The ecology and evolution of 
plant tolerance to herbivory. Trends in Ecology & Evolution 14:179— 
185. 

Thomas, H. Q., F. G. Zalom, and R. T. Roush 2010. Laboratory and field 
evidence of post-release changes to the ecological host range of Dio- 
rhabda elongata: has this improved biological control efficacy? Biolog- 
ical Control 53:353-359. 

Thornton, P. E., W. A. Kenney, and L. Zsuffa 1997. Generating surfaces 
of daily meteorology variables over large regions of complex terrain. 
Journal of Hydrology 190:214-251. 

Tracy, J. L., and T. O. Robbins 2009. Taxonomic revision and biogeogra- 
phy of the Tamarix- feeding Diorhabda elongata (Brulle, 1832) species 
group (Coleoptera: Chrysomelidae: Galerucinae: Galerucini) and 
analysis of their potential in biological control of tamarisk. Zootaxa 
2101:1-152. 

Turreson, G. 1930. The selective effect of climate upon the plant species. 
Hereditas 14:99-152. 

Vila, M., E. Weber, and C. M. Dantonio 2000. Conservation implications 
of invasion by plant hybridization. Biological Invasions 2:207—217. 

Weber, E., and B. Schmid 1998. Latitudinal population differentiation in 
two species of Solidago (Asteraceae) introduced into Europe. Ameri- 
can Journal of Botany 85:1110-1121. 

Whitham, T. G. 1989. Plant hybrid zones as sinks for pests. Science 
244:1490-1493. 

Whitney, K. D., R. A. Randell, and L. H. Rieseberg 2006. Adaptive intro- 
gression of herbivore resistance traits in the weedy sunflower Helian- 
thus annuus. American Naturalist 167:794-807. 

Wilcove, D. S., D. Rothstein, J. Dubow, A. Phillips, and E. Losos 1998. 
Quantifying threats to imperiled species in the United States. Biosci- 
ence 48:607-615. 

Williams, D. A., E. Muchugu, W. A. Overholt, and J. P. Cuda 2007. Col- 
onization patterns of the invasive Brazilian peppertree, Schinus tere- 
binthifolius, in Florida. Heredity 98:284-293. 

Williams, J. L., H. Auge, and J. L. Maron 2008. Different gardens, differ- 
ent results: native and introduced populations exhibit contrasting 
phenotypes across common gardens. Oecologia 157:239-248. 

Supporting Information 

Additional Supporting Information may be found in the online version 
of this article: 



392 



©201 4 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 381-393 



Williams et al. 



Invasive plant hybridization and biological control 



Appendix SI. ancova results for models investigating the effects of 
tamarisk introgression on various plant traits. 

Appendix S2. Tolerance scores and associated test statistics for each 
of the 43 plant genotypes involved in the outdoor garden study. 

Table SI. Site information, including latitude, elevation, frost-free 
days, annual extreme minimum temperature, and average species intro- 
gression levels for each of the tamarisk populations in the current study. 



Table S2. Competing models used in the Akaike's information crite- 
rion (AIC) model selection. 

Table S3. Results of Akaike's information criterion (AIC) model 
selection. 
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